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Abstract 

We study theoretically the behavior of a class of hydrodynamic dipoles. This study is motivated 
by recent experiments on synthetic and biological swimmers in microfluidic Hele-Shaw type geometries. 
Under such confinement, a swimmer’s hydrodynamic signature is that of a potential source dipole, and 
the long-range interactions among swimmers are obtained from the superposition of dipole singularities. 
Here, we recall the equations governing the positions and orientations of interacting asymmetric swimmers 
in doubly-periodic domains, and focus on the dynamics of swimmer pairs. We obtain two families of 
‘relative equilibria’-type solutions that correspond to pursuit and synchronization of the two swimmers, 
respectively. Interestingly, the pursuit mode is stable for large tail swimmers whereas the synchronization 
mode is stable for large head swimmers. These results have profound implications on the collective 
behavior reported in several recent studies on populations of confined microswimmers. 


1 Introduction 

Active systems, i.e., systems driven internally by self-propelled individual units, often exhibit rich collective 
behavior at the system’s scale; a scale that is typically several orders of magnitude larger than the scale of 
the individual unit. Such collective behavior naturally arises in disparate biological systems, from schools 
of fish [1] to suspensions of motile bacteria [2] and assemblages of sub-cellular extracts [3]. It also emerges 
in inanimate systems such as driven and self-propelled droplets and reactive colloids ma, and provide an 
attractive paradigm for reconfigurable smart materials and biomedical devices [?]• 

The question of how these highly-coordinated collective motions arise from piecewise interactions 
among individual units has been the subject of intense research in the past few years. A well-studied 
example is the behavior of self-propelled particles in a viscous fluid, [8]. Most of this work has focused on 
the instabilities and spatiotemporal fluctuations in three-dimensional (3D) systems. However, motivated 
by recent technological advances in producing and manipulating large ensembles of particles in microfluidic 
devices [SJIHlIin]; attention began to shift to the collective dynamics of particles confined in quasi two- 
dimensional (2D) geometries. Geometric confinement changes drastically the nature of the hydrodynamic 
interactions among particles, [9]. The long-ranged hydrodynamic interactions in 3D are driven by the force 
dipoles exerted by self-propelled particles on the fluid medium, [5]. In quasi-2D geometries, the solid walls 
screen the force dipole contribution, making it subdominant in comparison with the potential dipole arising 
from incompressibility, m- As a result, the long-range interactions among swimmers can be obtained from 
the superposition of dipole singularities. 

In this paper, we revisit the hydrodynamic dipole model proposed by Brotto et al. E] for asymmet¬ 
ric dumbbell swimmers in confined Hele-Shaw type geometries. The head-tail asymmetry causes a given 
swimmer to reorient, not only in response to the flow gradient as anticipated by Jeffery’s equation, but also 
in response to the flow velocity itself. This result is rooted in the fact that the lubrication forces between 
the swimmer and the solid walls hinder its advection by the fluid, inducing unequal translational motility 
coefficients at the swimmer’s head and tail. In m, Brotto et al. derived a kinetic theory-type model for 
a population of interacting swimmers and predicted a novel long-wave linear instability that leads to the 
emergence of large-scale directed motion and polarization in isotropic populations of confined large head 
swimmers. Lefauve and Saintillan m and Tsang and Kanso m used numerical simulations to explore 
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the implications of these instabilities on the collective behavior in finite-sized populations of interacting 
swimmers. 

The present paper examines the detailed dynamics of a pair of asymmetric (head-tail) swimmers in 
doubly-periodic domains, where the orientation dynamics is dominated by the flow field itself, thus neglecting 
reorientation in response to the flow gradient as done in mm- We particularly focus on a special class 
of solutions where the two swimmers move with constant speed and at constant orientation. We find two 
families of these “equilibrium-like” solutions: (1) both swimmers swim side by side in a synchronized way; 
and (2) one swimmer tailgates the other. We analyze their stabilities and find that they depend on the 
details of the head-tail asymmetry. We conclude this work by discussing the significance of these results to 
the behavior of populations of swimmers. 

Note that a dynamical theory of dipole interactions has also been pursued in two additional contexts. 
One motivation stems from the desire to obtain low-order representations of two-dimensional, inviscid and 
incompressible fluids in terms of interacting particles such as point vortices and point dipoles, see, e.g., [m- 
m- A shortcoming of these models is that the dipole’s self-propelled speed is ill-defined; thus, a dipole, 
unless properly desingularized, induces infinite velocity through its center. Another motivation for dipole 
models that is closer to the focus of this paper grew out of efforts to examine the role of hydrodynamic 
coupling in fish schooling. It is a well-known result in fluid dynamics that the leading order flow of a self- 
propelled body is that of a potential source dipole. Kanso and co-workers proposed a finite dipole dynamical 
system that captures the far-field hydrodynamic interactions of self-propelled bodies [EIEQ]. Each dipole 
consists of a pair of equal and opposite strength point vortices separated by a finite constant distance. By 
construction, the self-propelled speed is well-defined. These finite dipoles are advected by the local flow and 
reorient in response to the local flow gradient. In particular, the finite dipole reorients according to velocity 
gradient in the direction transverse to the dipole orientation, as opposed to reorienting in response to the 
flow gradient along the dipole’s direction predicted by Jeffery’s equation for slender bodies in viscous fluids. 
Intrigued by the similarities between the finite dipole model [Tll|20] and the dipole model of mm, Kanso 
and Tsang |22j presented a unified framework for deriving two point dipole models: a dipole consistent with 
the finite dipole model, appropriate for bluff bodies (fish) in potential flows, and another consistent with 
Jeffery’s equation for slender bodies and equivalent to the microswimmer model employed in mm- They 
further showed that, in unbounded domains, dipole pairs can synchronize their motions for a range of initial 
conditions; however, the details of the synchronized motion differ between the two models. 

The organization of this paper is as follows: In Section [21 we formulate the equations of motions for 
a systems of dipoles in unbounded and in doubly-periodic domains. A detailed treatment of the dynamics 
of dipole pairs is conducted in Section [31 These results are discussed in Section 0] in light of the large-scale 
simulations performed on the same system in [I21[T3] . 


2 Problem Formulation 


Microswimmer model. Consider a microswimmer composed of two connected disks of radii Rtaii and 
Rhead located at Ztaii and Zhead respectively, where z = x + iy is the complex coordinate (i = \/—1). Assume 
the two disks are connected by a frictionless rod of length i. The equations of motion for the swimmer’s tail 
{ztaii) and head (zhead) can be written in complex notation as (see [TTj l 


Ztail — U J- ytail'd^i^Ztail) j 

Zhead — UqB J- yhead’d^i^Zhead') ^headB 


( 1 ) 


Here, Uo is the swimmer’s self-propelled velocity, Oo its orientation angle, and w{z) is the velocity field of 
the ambient fluid. The bar notation denotes the complex conjugate, z = x — iy. The coefficients ytaii, y-head 
are the translational mobility coefficients whereas Xtaii, Xhead are unknown Lagrange multipliers that enforce 
the constraint \zhead — ztaii\ = In particular, the translational mobility coefficients ytaii and yuead arise 
from the balance of hydrodynamic drag and wall friction acting on the tail and head, and are decreasing 
functions of Rtaii and Rhead respectively, with values less than I, see, e.g. mm- 

We define the hydrodynamic center of the swimmer to be Zo = {XtaiiZhead + XheadZtaii)/{Xtaii + Xhead)- 
Our goal is to rewrite the equations of motion (HD in terms of the swimmer’s hydrodynamic center Zo and 
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orientation Oq. Let ^ ^ Rtaii, Rhead, and use Taylor series expansion to expand the flow velocity at the tail 
and head 


w{ztail) = w{zo) + 

'^{,^head) — “1“ 6 


, ^tail dw 

^tail ^head dz 
, ^head dw 

^tail ^head dz 


( 2 ) 


Substitute © into o to get that the equation governing the translational motion of the swimmer’s center 

Zo = Uoe~^°‘° + fj-w(zo), (3) 


where /r = {Xhead^J■tail + Xtail^J-head)/{Xhead + Xtaii)■ To obtain the equation governing the rotational motion 
of the swimmer, note that, by definition, £e'“° = Zhead — Ztaii, which gives, upon differentiating both sides 
with respect to time and further simplifications. 


do = Re 


{zhead 


ztajz)ie'“‘’ 

i 


(4) 


Here, Re denotes the real part of the expression in bracket. Now substitute m and into o to get 


do = Re 


dw. 

i/i—le 

dz 


2icto 




(5) 


where —— and w are evaluated at Zo and the constant parameters vi and 1^2 are given by 
dz 


Vl 


Xhead XfaH 


^2 — {fdhead fd,fail)/£. 


( 6 ) 


The sign of 122 dictates how the swimmer orients in local flow: it aligns to the local flow when U 2 > 0, 
that is, for large tail swimmers for which ^head — Idtaii > 0 (because fdhead/taii is a decreasing function of 
Rhead/taih [S])j and oppositc to the local flow when 1/2 < 0 , that is, for large head swimmers for which 

f^head Idfail ^ 0 - 


Hydrodynamic interactions of multiple microswimmers. Consider the interaction of multiple mi¬ 
croswimmers in an unbounded fluid domain. By virtue of ([3]) and ®, the dynamics of N swimmers, all 
having the same self-propelled velocity U, can be expressed in concise complex notation 


Zn = Ue -I- fiw{zn), 


cuyj — R,e 


dw. 

1^1—le 
dz 


2ictn 


+ i22wie^°‘ 


(7) 


Here, Zn and denote the position and orientation of each swimmer (n = 1,..., N). To close the model, 
one needs to obtain an expression for the fluid velocity field w{z). Recalling that each swimmer induces 
a far-held velocity which is that of a potential source dipole [9], the far-held how of a microswimmer j 
located at Zj = Xj + \yj and oriented at an arbitrary angle Uj can be described by the complex velocity 
w{z) = Ux — iuy = l{z — Zj)'^, where a is the dipole strength. Note that a = R^U, where R is the 

effective radius of the swimmer. A microswimmer n responds to the how induced by all microswimmers in 
the huid domain, namely, 


N 

w{Zn) = X] 

J = 1 


e J 


{Zn Zj)"^ 


( 8 ) 
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Figure 1: (a) Two modes of solutions are obtained from llllll ; a pursuit mode (blue) where the swimmers trail one another 

and a synchronization mode (red) where they swim side by side. The shown solid curves correspond to coi = —iuj 2 = 5 and the 
dashed lines to lji = —1012 = 10. The separation distance c is set to c = A. (b) Summary of the stability analysis for these two 
modes. 


Microswimmers in doubly-periodic domains. When the swimmers are placed in a doubly-periodic 
domain, one needs to take into account, not only the velocity field induced by the swimmers themselves 
but also the effect of their image system. A given swimmer n has a doubly-infinite set of images. Thus, 
evaluating w{z) requires the evaluation of conditionally-convergent, doubly-infinite sums of terms that decay 
as l/l^p. These sums are evaluated using an approximate numerical approach in [12) . In [13) . we offered a 
closed-form analytic expression for these infinite sums in terms of the Weierstrass elliptic function, namely, 

N 

w(z) = ^ crp(2: - Zn; Wi, W2)e‘“". (9) 

n—1 

The Weierstrass elliptic function p[z) is given by p(z;wi,W 2 ) = ^ + Thki ( (z-hki)'^ ~ Vtki = 

2k(jj\ -\- 21 uj 2, k,l £ Z—{0}, and uji and UJ 2 being the half-periods of the doubly-periodic domain. This 
function has infinite numbers of double pole singularities located at 2 ; = 0 and 2 : = flfc;, corresponding to 
the l/\z\^ singularities induced by the potential dipoles. Equations ([7]) and (O form a closed system for N 
swimmers in a doubly-periodic domain. 

We conclude by writing the system of equations © and (|ni) in dimensionless form using the swimmers 
radius R as a length scale and R/U as a time scale. That is, we introduce the dimensionless spatial variable 
z = z/R and time variable t = tU/R. We then drop the tilde notation assuming all variables are non- 
dimensional. Equations o and m have the same form but the parameters U and a are now normalized to 
one, that is, t/ = 1 and (7 = 1. The parameter values p, vi and 1^2 are also non-dimensional. 


3 Pursuit and synchronization 

We consider two microswimmers in a doubly-periodic domain, and focus on their dynamic response when 
vi = 0, that is, when their alignment with the flow gradient is negligible. In this case, the orientation 
dynamics is dominated by alignment with the flow due to head-tail hydrodynamic asymmetry. 

Periodic solutions and relative equilibria. We look for special solutions where the two swimmers move 
at the same velocity and orientation for all time. That is, we look for solutions where zi = Z 2 = constant 
and di = 02 = 0. To obtain the initial conditions that lead to this behavior, it is convenient to rewrite the 
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Figure 2: Aperiodic behavior of two dipoles in doubly-periodic domain. The parameter values are a = tt/S, -Zi(O) = —2exp{i6), 
22 ( 0 ) = 2exp{i9), while 0 is obtained by numerically solving the first condition in l IllD . The positions of the dipoles are marked 
by ‘x ’ at t = 0 and by ‘o’ at the end of the integration. As time progresses, the two trajectories densely fill the whole domain. 


equations of motion dZH) in terms of the reduced coordinate zi — Z 2 which we set to zi — Z 2 = /3 = ce*® (see 
inset of Figure (Ha)). To this end, one gets 


/3 = e"'“^ - e"‘“= + ^p(/3)(e*“= - e*“^), 

di = d2 = J^2R.e[ie*^“^+“=V(/3)]- 


( 10 ) 


The translation equation for P is identically zero when ai = a 2 = a. Whereas to guarantee di = d 2 = 0, 
one must satisfy the condition 


Re[/ci(ce ‘®)] = 0 therefore “ “ ^ (o^i 

or 


Im[yo(ce '®)] 
. Re[p(ce“*®)] 


tan2Q;, Re[p(ce '®)] 0, 


(a, a) 

(a, a + 7r/2) 


( 11 ) 


A total of ten strict relative equilibria of the two swimmers are depicted schematically in Figurelija). Namely, 

TT TT StT 

the five solutions given by a = 0, —, —, — , 7 r and 9 = a correspond to the two dipoles moving parallel to 

TT TT StT 

each other in a “pursuit” mode (blue arrows), whereas the five solutions given by a = 1 *’ "J’ 

0 = a + 7r/2 correspond to the two dipoles moving in tandem in a “synchronized” mode (red arrows). 

The second set of solutions, that is, the values of (a, 0) for which Im[p(ce“‘®)]/Re[yo(ce“'®)] = — tan2Q; 
and Re[p(ce“'®)] 7 ^ 0, are not analytically available and need to be computed numerically. Figure [lla) 
shows the values of (a, 0 ) that satisfy these conditions - clearly, two branches of solutions are obtained. 
These solutions depend implicitly on the domain size (wi, ^ 2 ) and on c, the separation distance between the 
two swimmers. In other words, for a choice of domain size and separation distance c, (a, 0) are computed 
accordingly such that the two dipoles move at the same constant velocity and orientation for all time. The 
two branches shown in Figure[T](a) correspond to two modes of behavior: a pursuit mode where one swimmer 
trails the other, and a synchronization mode where the two swimmers move side by side. These solutions, 
while they correspond to the dipoles moving at constant velocity and orientation, can exhibit two distinct 
types of dynamical behavior due to the doubly-periodic nature of the domain, namely, they could lead to 
aperiodic and periodic motion of the dipoles. Aperiodic motion refers to the case where the paths of the 
dipoles densely fill the whole domain, as shown in Figure [2] This seems to be the generic behavior for 
arbitrary initial conditions. Periodic behavior refers to trajectories that satisfy the condition 


zi(r) = 2 : 1 ( 0 ) -I- 2pa;i -b 2ga;2, Z2(T) = 2 : 2 ( 0 ) + 2puji + 2qu}2 


( 12 ) 
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Figure 3: (a) Pursuit and (b) synchronization in two dipoles undergoing periodic motion. The parameter values are a = 
tan“^(3), zi{0) = —2ex.p{i6), 22 ( 0 ) = 2ex.p{i0), while 6 is obtained by numerically solving the first condition in JllH . The 
positions of the dipoles are marked by ‘x ’ at t = 0 and by ‘o’ at the end of the integration time t = 80. 




Figure 4: Pursuit and synchronization modes as attracting modes. Two dipoles hone in on quasi periodic trajectories where: 
(a) one dipole is in pursuit of the other for U 2 = 0.5. (b) the two dipoles synchronize and move along parallel trajectories for 
V 2 = —0.5. The initial conditions are ^i(O) = Q:i(0) = 0 : 2 ( 0 ) = 0 while 22 ( 0 ) = 1.5 + li in (a) and 22 ( 0 ) = 2 + li in (b). 


where p and q are integers and T is the period of the motion. This amount to the additional condition 

ai(0) = 02 ( 0 ) = tan“^(-). (13) 

P 

The ratio oi q/p indicates the ratio of the number of times the dipole crosses the y and x axes in one period 
T. Figure [3] depicts the periodic behavior of two dipoles in pursuit and synchronization modes for q/p = 3. 


Stability analysis. We analyze the linear stability of the pursuit and synchronization modes by considering 
small perturbations <5/3 = 6Px + i^/3y, 5oii and 5a2 about /3 = ce‘® {Px = ccos0, and Py = csin0) and 
«! = 02 = ot, with (a, 6) satisfying (fTTl) . We linearize equations (flUl) accordingly. The linearized equations 
can be written in matrix form as follows: 



( SPx \ 


( SPx \ 

d 

SPy 

= M 

SPy 

dt 

Sai 


5ai 


\ (502 ) 


\ (50(2 / 
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where the Jacobian matrix M is given by 


( ^ 

0 

:/2Re[ie‘^“p'(/3)] 
y :/2Re[ie‘^“p'(/3)] 


0 

0 

-j^2R.e[e‘^“p'(/3)] 

-j^2R.e[e‘^“p'(/3)] 


— sin a — /rRe[ie*“p(/3)] 

— cos a — ^Im[ie‘“p(/?)] 

-:/2Re[e*^“p(/3)] 

-:/2Re[e*^“p(/9)] 


sina + ^Re[ie‘“p(/3)] \ 
cos a + /rIm[ie*“p(/3)] 
-^2Re[e*^“p(/3)] 
-^2Re[e*^“p(/3)] / 


(15) 


We compute the eigenvalues numerically and find that, for large tail swimmers > 0, the pursuit mode 
is stable, whereas for large head swimmers V 2 < 0, the synchronization mode is stable. Our findings are 
summarized in Figure [TJb). 

We test our results numerically by integrating the nonlinear equations da© for arbitrary choices 
of initial conditions. Interestingly, the pursuit and synchronization modes seem to be globally attracting 
modes in the case of large tail and large head swimmers, respectively. Figure ©a) shows a depiction of 
two large-tail swimmers honing in on quasi-periodic pursuit trajectories, while (b) depicts two large-head 
swimmers synchronizing their motion in finite time to swim side by side. 


The limit of unbounded domain. We conclude this section by noting that in the limit of infinite domain, 
the solutions m of the doubly-periodic system da© converge to the relative equilibria of the unbounded 
system da©- In the unbounded system, the relative equilibria can be obtained either by symmetry arguments 
or by analytical manipulation of the equations of motion. Namely, one has two families of relative equilibria 
9 = a and 9 = a -\- ii 12 which correspond to pursuit and synchronization trajectories, respectively. The 
convergence of the solutions in (ED to these solutions is relatively fast, as indicated in Figure ©a). As wi, 
a ;2 —>■ oo, the Jacobian matrix M converges to 


= 


2v2 

c? 

2v2 


0 

0 

sin(2a — 39) 


2v2 
2v2 


0 

0 

cos(2a — 39) 


, sin(2a — 30) cos(2a — 30) 

\ C'’ 


— sin a -I- sin(a — 20) sin a — sin(a — 20) \ 

— cos a — ^ cos(q; — 20) cos a^ cos(a — 20) 

— ^ cos(2q! — 20) — cos(2a — 20) 

c‘^ &■ 

— ^ cos(2q! — 20) — ^cos(2a —20) . 


(16) 


The corresponding eigenvalues are [0, 0, 0, =f2i^2/c^]- The eigenvalue —2v2l(? corresponds to the pursuit 
mode where a = 9, whereas -\-2v2l(? corresponds to the synchronization mode. This means that, for large 
tail swimmers with U 2 > 0, the pursuit mode is linearly stable and the synchronization mode is unstable, 
whereas for large head swimmers when V 2 > 0, the opposite is true, thus confirming the results obtained 
above for finite-sized doubly-periodic domains. 


4 Discussion 

We revisited the hydrodynamic dipole model governing the interaction of asymmetric microswimmers in 
Hele-Shaw confinement, El- Following [T3l[20], we obtained a closed-form expression for the velocity field 
induced by the swimmers and their image system in doubly-periodic domains. We treated in details the 
dynamics of two interacting swimmers, and found two special solutions that correspond to pursuit and 
synchronization of the two dipoles. The pursuit mode is stable and attracting for large tail swimmers while 
the synchronization mode is stable and attracting for large head swimmers. By attracting, we mean that, 
starting from arbitrary initial conditions, large tail swimmers tend to tailgate each other while large head 
swimmers tend to synchronize their motion in finite time to swim side by side. These results are particularly 
interesting in light of the collective behavior reported in [I21I13] on populations of such swimmers. In these 
works, large tail swimmers were observed to “develop active lanes” E] “tail-gate each other” [13], as 
shown in Figure©top row), which suggests that the pursuit mode remains stable as the system size increases. 
Note that, to generate Figure© we integrate equations (© and (© for a population of 400 dipoles, starting 
from a uniform isotropic distribution and using the parameter values /i = 0.9, = 0 and i /2 = 1 (large tail) 

or V 2 = —0.5 (large head). 
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Figure 5: Emergent collective behavior in large tail swimmers (top row) and large head swimmers (bottom row) starting from 
a uniform isotropic distribution. Large tail swimmers tend to tailgate each other, thus forming active lanes, while large head 
swimmers tend to form stationary clusters. Parameter values are p = 0.9, z /2 = 1 (top row) and V 2 = —0.5 (bottom row) in a 
total population of 400 swimmers. 


Populations of large head swimmers were shown to form heavily polarized sharp density waves in m, 
consistent with predictions based on linear stability analysis of a kinetic-type continuum model m- One 
could conjecture that the synchronization mode observed here in pairs of large head swimmers may be 
responsible for the global polarization observed in m- However, this thinking is too simplistic. The 
emergence of global polarization patterns in finite size populations is not intuitive given the nature of dipolar 
interactions among the swimmers. Further, these polarized density waves were not observed in the detailed 
parametric study reported in m Figure 7]. Instead, [13] reported, in agreement with unpublished results 
by Levaufe and Saintillan, that large head swimmers tend to form stationary clusters (see Figure [5l(bottom 
row)), which are not predicted by the linear stability analysis of [IT]. All this is to say that the global patterns 
of the finite size systems in |12j and |13| are in agreement, except for the global polarization pattern. This 
inconsistency may be due to differences in the system size - thousands of particles in |12| versus hundreds in 
Figure [5] and in m - or to differences in the details of the numerical implementation. In |12] , the point dipole 
model is desingularized and hydrodynamic interactions are approximated for fast computations, whereas m 
use a local repulsion potential for collision avoidance and accurately account for hydrodynamic interactions 
and the doubly-infinite image system. While the difference in numerical implementation may play a role, the 
system size may be the main reason why polarized waves are not observed in m- Brotto et al. m predicted 
this behavior for a continuous kinetic-like model, therefore it is not surprising that it is not reproduced by 
a fully nonlinear model with only a few hundred swimmers. Irrespective of the reason, the results reported 
in this study suggest that the global polarization mode in large head swimmers is not “robust” to system 
perturbances, whereas the pursuit mode in large tail swimmers is. 
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